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Abstract. We present the first attempt to determine the iron 
abundance in hot central stars of planetary nebulae. We perform 
an analysis with fully metal-line blanketed NLTE model atmo- 
spheres for a sample of ten stars (T e gp ~ 70 000 K) for which 
high-resolution UV spectra are available from the IUE archive. 
In all cases lines of Fe VI or Fe VII can be identified. As a gen- 
eral trend, the iron abundance appears to be subsolar by 0.5- 
1 dex, however, the S/N of the IUE spectra is not sufficient 
to exclude a solar abundance in any specific case. Improved 
spectroscopy by either FUSE or HST is necessary to verify the 
possibility of a general iron deficiency in central stars. The sus- 
pected deficiency may be the result of gravitational settling in 
the case of three high-gravity objects. For the other stars with 
low gravity and high luminosity dust fractionation during the 
previous AGB phase is a conceivable origin. 

Key words: stars: abundances - stars: atmospheres - stars: 
evolution - stars: AGB and post- AGB - white dwarfs - Ultra- 
violet: stars 



1. Introduction 

The continuing spectroscopic study of central stars of plane- 
tary nebulae (CSPN) is motivated by two observational facts 
which reveal our incomplete understanding of post- AGB stellar 
evolution. The first concerns the mere existence of hydrogen- 
deficient objects which comprise the spectroscopic types [WC] 
and PG1159. The second concerns the apparent lack of 
hydrogen-rich CSPN at the hot end of the white dwarf cooling 
sequence, suggesting that during this stage all H-rich objects 
become H-deficient. The first problem initialized an almost 
complete analysis of all known H-deficient post- AGB stars, re- 
sulting in accurate determinations of the effective temperature 
(T c ff) and surface gravity (g), as well as the abundance of the 
main atmospheric constituents (He, C, O) and in some cases 
trace of elements (H, N, Ne). The second problem could be 
solved by the systematic search for (pre-) white dwarf central 
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Fig. 1. Location of the program stars in the logT c ff - 
\ogg diagram compared with post- AGB evolutionary tracks 
(jSchonbei-ner 1983|; [Blocker & Schonberner 1990| ) for H- 



burning stars with different remnant masses (in Mq) as indi- 
cated. 



stars. Subsequent analyses of the H-rich objects revealed T e g, 
log g, and H:He ratio and the results complement the analyses 
of other H-rich CSPN performed in numerous previous analy- 
ses (see Tab. ||). 
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Table 1. The program stars and their atmospheric parameters. Values for T c g [K], log g [cm s~ 2 ], and helium abundance were 
taken from literature. Columns 5 and 6 list our results for the Fe abundance and the H I column density along the line of sight as 
derived from the Lya profile; the references are listed in Tab. [2] 



PN 




logs 


nHe/nHeQ 


log nFc/nFcQ 


nm 


[cm- 2 ] 


references 




SWP numbers of co-added spectra 


Abell 36 


93 000 


5.3 


1.5 


— 1.3 


1.9 


■10 19 


4, 5, 8, 9 




16478, 


38132, 


47814, 47815 


LSE 125 


85 000 


5.1 


0.5 


— 0.5 


3.0 


■ 10 19 


2,5 




30269 






LSS 1362 


100000 


5.5 


1.0 


-0.5 


1.6 


■ 10 19 


3,4,5,8 




19830 






NGC 1360 


110000 


6.0 


1.0 


-0.5 


3.0 


■ 10 18 


1,2, 4,5,* 


», 10 


55902, 


56037, 


56038 


NGC 1535 


70000 


4.6 


1.2 


-0.5 


1.4 


■10 19 


1,2, 4,5,* 


! 


56065, 


56066 




NGC 4361 


82 000 


5.5 


0.5 


-1.3 


5.0 


■10 18 


1, 2, 4, 5, i 


,.9 


20440 






NGC 6853 


99 000 


6.8 


1.0 


-0.5 


3.0 


■10 18 


5, 7, 8 




18340 






NGC 7009 


82 000 


4.8 


0.5 


-0.5 


1.5 


■ 10 19 


2, 4, 5, 8 




05174, 


08580, 


08680, 23383 


NGC 7293 


107 000 


7.0 


0.5 


-0.5 


5.5 


■10 18 


4,9 




19859, 


48055, 


48063 


S216 


85 000 


6.9 


0.1 


-0.5 


5.0 


■10 18 


4, 6, 7, 8 




56071, 


56072 





Table 2. References in Tab. [j] 



5F 



Vlendez et al. (1981| ) 
Heberetal. (1988|) 
Vlendez (199 1|) 
^apiwotzki (19931) 
juigley et al. (1993[) 



Mendez etal. (1988) 



Feibelman & Bruhweiler (1990b 



Tweedy & Napiwotzki (1992; ) 



8 Tweedy (1993]) 



10 Hoareetal. (1996) 



Almost no effort has been yet made to analyze metal abun- 
dances in the H-rich CSPN. Interesting insight into AGB evo- 
lutionary phases may be gained by such analyses, because ob- 
vious abundance anomalies exist among stars w ith otherwise 
similar atmospheric parameters ( Mendez 1991 ). Before em- 
barking on such a project and discussing the abundance pat- 
terns it is reasonable to obtain information about the primordial 
metallicity of these objects by analyzing their iron abundance. 
The reason herefore is that iron is not affected by nuclear pro- 
cesses during previous evolutionary phases. Such an analysis 
requires UV spectroscopy. As a first step it is mandatory to 
check, if the IUE Final Archive contains useful data. It is the 
aim of this paper to exploit that archive exhaustively in order 
to identify and analyze quantitatively iron lines in hot central 
stars. 



2. Sample selection and IUE data reduction 

Our sample selection started from a CSPN compilation by 
Mendez (1991 ) and we complemented his list of H-rich objects 
by several objects dealt with in the more recent literature (e.g. 
see Tab. ^|). For these stars we checked the IUE archive for high 
resolution SWP spectra and found data for 27 H-rich CSPN. 
From this sample those objects were excluded, which have high 
mass-loss rates, i.e., those showing Of-type optical spectra and 
those with extraordinarily strong P Cygni resonance line pro- 
files in the UV. Thereby the final sample consisted of ten stars, 
which can be analyzed reliably with our static model atmo- 
spheres. All program stars have T c g ~ 70 000 K. Fig. filshows 
their position in the log T ff — log ^-diagram and Tab!"]]] sum- 
marizes their atmospheric parameters. 
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Fig. 2. Calculated continuum of a typical co-added IUE spec- 
trum 

The SWP spectra were extracted from the IUE Final 
Archive and subject to several reduction steps before compari- 
son with model spectra. These steps comprise the assembly of 
single Echelle orders in order to obtain a complete spectrum 
in the 1200- 1600 A range; elimination of the strongest inter- 
stellar absorption lines by Voigt profile fitting (Lya and Si II 
lines); correction for radial velocity; normalization of the spec- 
tra. The last step is rather delicate. Firstly the ripple correction 
for the Echelle orders in the archival data is far from being per- 
fect (see Fig. |2|) and has to be corrected. Secondly, fixing the 
continuum is a rather subjective procedure and thus the major 
error source for the quantitative analysis, because of the gen- 
erally high noise level in the spectra. To determine the con- 
tinuum we smoothed the spectrum with a median filter of 2 A 
width, followed by a convolution with a Gaussian profile of 2 A 
FWHM. This method is inapplicable in regions with P Cygni 
profiles and the continuum has to be determined manually. Fur- 
thermore the determined continuum around deep ISM lines is 
too low and has to be raised manually. 
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Fig. 3. Distribution of the analyzed CSPN (O) and models (□) 
in the T c g - log g - plane. For each point of the grid spectra for 
eight different Fe and Ni abundances have been calculated 
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Fig. 4. Ionization structure of iron in a model atmosphere with 
T off = 110000 K, logg = 6.0 and solar Fe abundance. Note 
particularly the relative occupation of Fe VII and Fe VIII in the 
deeper layers of the atmosphere where the iron lines are formed 
(logm = -2...0) 



3. Model Atmospheres 

We have calculated a grid of plane-parallel non-LTE model 
atmospheres (Fig. ^ in radiative and hydrostatic equilibrium. 
The computer code is based on the Accelerated Lambda Itera- 
tion method (Werner & Husfeld 1985; Werner 1986) and it can 



handle the line blanketing of iron group elements by a statisti- 
cal approach using superlevels and superlines with an opacity 



sampling technique (Anderson 1985; Dreizler & Werner 1993) 



The latest version of the code and a detailed description of the 



numerical method can be found in Werner & Dreizler (1999). 
Our model atmospheres include the CNO elements as well as 
iron and nickel self-consistently, i.e. their back-reaction on the 
atmospheric structure is accounted for. In es sence, they are ver y 
similar to the models described in detail by fr-taas et al. (1996 ); 



Table 3. Summary of model atoms used in our model atmo- 
sphere calculations. Numbers in brackets denote individual lev- 
els and lines used in the statistical NLTE line blanketing ap- 
proach for iron and nickel. The model atom for each chemical 
element is closed by a single level representing the highest ion- 
ization stage (not listed explicitly) 



element 


ion 


NLTE levels 


lines 




H 


I 


1 ft 




36 




He 


I 


5 




3 






II 


10 




36 






III 


6 




A 






T V 
1 V 


4 




1 
1 




N 


IV 


6 




4 






V 


4 




1 




O 


IV 


1 











V 


6 




4 






VI 


4 




1 




Fe 


IV 


7 


(6472) 


25 


(1027 793) 




V 


7 


(6 179) 


25 


(793 718) 




VI 


8 


(3 137) 


33 


(340 132) 




VII 


9 


(1 195) 


39 


(86 504) 




VIII 


7 


(310) 


27 


(8 724) 


Ni 


IV 


7 


(5 514) 


25 


(949 506) 




V 


7 


(5 960) 


22 


(1006189) 




VI 


7 


(9 988) 


22 


(1 110 584) 




VII 


7 


(6 686) 


18 


(688 355) 




VIII 


6 


(3 600) 


20 


(553 549) 


total 




135 


(49041) 


346 


(6565 054) 



Haas (1997 ), so we restrict ourselves here to a summary of 
the model atoms in Tab. || The main extension of that mod- 
els refers to a detailed model atom for Fe vm. This turned out 
to be essential in order to reliably compute Fe VII lines, because 
the relative occupation of both ionization stages sensitively de- 
pends on the detailed treatment of their respective level popu- 
lations via ionization and recombination processes. 

As an example, the iron ionization structure in a typical 
model atmosphere is shown in Fig. Q 



4. Results and discussion 

We derived our results with two independent methods. First we 
compared the observed and theoretical spectra by eye (Fig. |]), 
second a x 2 -test has been performed (Fig. ||). To get more re- 
liable results, the x^ ci has been computed only around strong 
iron lines excluding known ISM lines. As a general trend, the 
iron abundance appears to be subsolar by 0.5-1 dex with both 
methods (see Tab. [j]). However, the S/N of the IUE spectra is 
not sufficient to exclude a solar abundance in any specific case. 

We cannot identify nickel lines for certain, which corre- 
spond to an upper limit of a solar nickel abundance. Ni V lines 
have been detected in other central stars and sdO stars by 




Fig. 5. Details from fits to the iron line spectrum of two central stars. Top: S 216; T cS = 90 000 K, log g = 7.0; Bottom: NGC 1360; 
T c ff = 1 10 000 K, log g = 6.0. Models (dotted): a solar Fe and Ni; b 0.5 dex subsolar Fe and Ni; c comparison of both models; d 
strong Fe VI and Fe VII lines are marked. 
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Fig. 6. Contour plot of the \ 2 test for the spectrum of NGC 1360. In the left plot the x 2 cd values for those models with logg at 
6.0 and T c g around 110 000 K (40 = 8-5 models, see Fig. |]) are shown, whereas in the right plot T e g remains fixed and log g 
varies. The high \ 2 values (> 5) are due to the low S/N level and the imperfect normalization (see Fig. ||). As a general trend, 
the x 2 test favors models with subsolar iron abundance 



Schonberner & Drilling (1985), however, these stars have ef- 
fective temperatures which are lower than those of our objects. 

For three high-gravity objects in the sample the suspected 
iron deficiency may be the result of gravitational settling. For 
the other stars of the sample with low gravity and high lu- 
minosity another explanation must be found since mass loss 
efficiently suppresses gravitational settling (Unglaub & Bues 
1998). Following suggestions in the literature (Lambert et al. 
1988; |Venn & Lambert 199$ ftond 1991| ; Van Winckel et al. 
1992), the iron abundance in these stars might not be primor- 
dial, but the result of chemical fractionation through dust for- 
mation and subsequent loss due to stellar wind during previous 
AGB phases. 

Therefore improved spectroscopy either by FUSE or HST 
is necessary to confirm or reject the possibility of a general Fe 
deficiency in these central stars. Furthermore, the determina- 
tion of the abundance of CNO-elements is desirable in order to 
get a better understanding of the chemical evolution of these 
stars. 

Our study suggests that iron is not necessarily a suitable 
tracer to determine the primordial metallicity. Alternatively one 
may use the S or Zn abundance to obtain reliable information 
(Van Winckel et al. 1992). The abundance of both elements is 



not affected by nuclear processes during previous evolutionary 
phases and both elements do not tend to dust formation as iron 



does (Habing 1996 ). Due to the lack of detailed atomic data for 
Zn, one has to concentrate on the S lines in the FUV and UV 
region. This emphasizes the need for improved spectroscopic 
data to solve these questions. 
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